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I. INTRODUCTION 



The purpose of this paper is to examine the feasibility 
of using a micro computer in Naval aircraft to aid the pilot 
in the solution of the ballistics problem for unguided 
weapons. Specifically, this problem is to release a weapon 
at some point in space with the proper velocity vector so 
that the weapon will eventually strike a selected target. 
Since several models of attack aircraft (A-6E, A-7E} are 
equipped with computers which perform this function, special 
attention will be focused on those aircraft in which the 
pilot manually provides the solution. Such aircraft are 
referred to as fixed sight aircraft and the F-4 is 
representative of this class. 



The micro computer, because of its low cost an 
is a device with the potential to provide the pilot 
sight -aircraft with useful and timely informa 
determine whether such a potential could be exploi 
INTEL 8080 is examined as being representative 
micro computers and the F^-4 is considered as a typi 
sight aircraft. The concepts, however, are in ten 
more general in nature and no discussion , 
conclusion contained herein is intended as an endor 
criticism of the manufacturer or weapon system invo 
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II 



UKGUIDED 



WEAPONS DELIVERY USING FIXED SIGHT AIRCRAFT 



A. DEFINITIONS 



Fixe 
equipped 
computer 
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an angle 
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range of 
they rem 
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n below 
REFERENCE 
Typically 
zero to a 
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rcraft is a term applied tc those 
ther a lead computing gunsigh 
bomb release device. In such a 
r gunsight is located in the 
iew. Its purpose is to represent 
adians (mils) measured in the 
a longitudinal reference line c 
LINE (FEL) . Figure 1 depicts 
these devices are adjustable 
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til manually changed by the pilot 
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Fig 1 SIDE VIEW AND PILOT VIEW OF A GUNSIGHT 



10 



B • 



PROBLEM DEFINITION 



For a particle of constant mass, m, moving under the 
action of an applied force, ?, Newton’s second law of motion 

states that 



F=m*a. {1} 

If the weapon is considered to be such a particle, and 
if the position of the weapon is measured in a rectangular 
coordinate system, (x,y,z), then equation (1) may be written 
in vector form as 

(a ,a , a )= (F ,F ,F )/m. (2) 

x y z x y z 



For air to ground 
system is chosen so that 
the down range travel of 
displacement in azimuth, 
the moment of release, x 



weapons delivery, the coordinate 
y represents altitude, x represents 
the weapon, and the z axis measures 
The origin is chosen so that, at 
equals zero. 



If it is assumed that the air mass through which the 

weapon is moving is stationary with respect to the earth, 

that gravitational forces are constant with respect to y, 

and that forces induced by the rotation of the earth are 

negligible, then F =0. The remaining forces acting on the 

z 

weapon are drag, thrust, and gravitational force. It will 

be assumed that thrust is zero, for all weapons except 
rockets. Rockets are assumed to have a short initial period 
of constant thrust after which they fall with the same 
discipline as all other unguided weapons. Drag acts in the 



1 1 



direction opposite to the direction of the velocity vector 

2 

and is directly proportional to !|V|| and air density. In 

a standard air mass, density is a decreasing function of 

altitude. Therefore equation (2) may be written in 
component form as follows 



a =-drag * 1/m (3) 

x x 



a = -(drag +g*m)/m (4) 

Y y 



where g is a constant 32.17 ft/sec . 



Letting THETA=ARCTAN (v /v ), (3) and (4) become 

y x 

a =-drag*cos (THETA) /m (5) 
x 

and 



a =-drag*sin (THETA) /m -g. (6) 

y 



Drag is represented as 



drag=rho*S*V *C 



where * rho* is the air density, ’S' is the cross sectional 



2 2 1/2 

area, ' V* is the total velocity ( (V +V ) 

x y 



) , and C , 

d 



the drag coefficient, is a dimensionless factor whose value 



depends on the shape of the weapon and the ratio of *V' to 
the speed of sound at altitude y. Equations (5) and (6) may 
be considered as differential equations 
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a =f (v , v , x , y, t) 
x x y 



(7) and 



a =h (v 

y x 



/V ,x,y,t) 

y 



( 8 ) 



In air-to-ground weapons delivery, equations (7) and (8) 
may be solved numerically for any set of initial conditions 
0 0 0 0 

(v ,v ,x ,y ) to determine a time t' where y =the target 
x y t * 



elevation and x 

t ' 



is the down range travel of the weapon. 



Therefore, the problem for air-to-ground delivery is to 
position the aircraft so that the initial conditions are 

matched when the selected aim point is a distance x down 

t ' 



range from the aircraft. 



For unguided air-to-air weapons, equations (7) and (8) 
may be solved together with initial value conditions to 
describe the position (x (t) , y (t) , z (t) ) of the weapon at any 
time t after the weapon has been fired. Since the position 
of the target is also some function of time, 

(x (t),y (t),z (t) ) , the problem is to find a 

target target target 

set of initial value conditions such that for some time, 

f ' , 

<x(f) ,y(t”) ,z(f)) = (x (t * • ) , y (t * * ) , z ( 

target target target 

t")) • 



In either case, the aircraft instruments provide the 
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dive 



initial velocity vecto 
information is not directl 
instrument and the inform 
(true airspeed and altitude 
the pilot before he can 
further delay. Finally, th 
instantaneously to control 
that sudden perturbations o 
pilot induced changes in p 
will cause the information 
be unreliable for several 
in Fig 2, the gunsight ima 
the earth so that the 
precisely once during the d 



r . Furthermore, dive a ngle 
y available from any aircraft 
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) must be mentally processed by 
physically react , thereby adding 
e aircraft itself cannot react 
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f the angle of attack caused by 
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provided by a fixed gunsight to 
seconds. In addition , as shown 
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proper sight 'picture' occurs 
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Fig 2 SIGHT LINE MOTION IN A DIVE 
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In spite of these difficulties, there are strategies 
which help the pilot to learn a solution to the problem. 
First, the deliveries have been standardized so that the 
pilot may concentrate his efforts at mastering several basic 
maneuvers. This practice tends to reduce the magnitude of 
delivery errors. The second aspect which aids the pilot is 
known as error analysis. Since actual conditions rarely 
match those cf the practice patterns, the pilot is taught to 
recognize those errors which are present and to apply 
corrective measures. Fig 3 graphically depicts the effect 
of various errors in airspeed, dive angle, release altitude, 
and aim point. For each of the standard deliveries, there 
is a table in reference which lists, for each weapon, 
the effect upon weapon impact of errors of 10 knots in true 
airspeed ( fast and slow ) , one degree in dive angle ( steep 
and shallow) , 100 feet in release altitude ( high and low) , 
and one raiiliradian in sight setting ( plus and minus) . 



To use the tables, the pilot memorizes key v 
the weapons he is going to release. During th 
portion of his dive, ne first attempts to min 

errors in each of the delivery parameters. If 

that the value of some or all of the delivery 
will not equal the preselected valuer , he est 
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causes another error either in release c ltitude 
angle to counteract the effect of the previous err 
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As an illustration of this method, suppose the pilot 
estimates that he will be 10 knots fast and one degree steep 
at release altitude and that the combined effect of the 
errors will place the weapon 100 feet past the target. If, 
in the error sensitivities chart for the particular weapon, 
the entries for +100 ft. and -1 mil are -25 feet and -10 
feet respectively, then the pilot may choose to release the 
weapon 400 feet high in relation to the preselected release 
altitude or he may use an imaginary sight line which 
represents an angle of 10 mils less than the preselected 
angle . 

In view of the difficulties involved in fixed sight dive 
bombing, it is not surprising that considerable effort has 
been devoted to the development and implementation of 
computer systems to improve the delivery accuracy of all 
current generation Naval Attack aircraft. 



D. AIR TO AIR GUNNERY 

In air-to-air gunnery, the pilot of a fixed sight 
aircraft essentially solves his aiming problem in the same 
manner as for air to ground weapons delivery. He selects 
several representative situations and consults a table to 
determine the correct sight setting for each particular 
situation. In air to air gunnery, however, the pilot would 
usually like to be able to fire his aircraft's guns at any 
time he has a good chance of hitting the target. Therefore, 
rather than try to match any preselected set of firing 
conditions, the pilot takes the best guess and fires. The 
guns are loaded with a mixture of tracer rounds so that the 
pilot visually adjusts his aim depending on his evaluation 
of the effect of previous aim angles. 
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PROPOSED SOLUTION 



In order to achieve the full bene 
the micro computer, it was assumed tha 
equipment such as doppler, air-to-grou 
inertial navigation devices would 
modified system. Instead, the sys 
around existing components with ultima 
problem solution still resting with 
computer is used chiefly to perform th 
the pilot of fixed sight aircraft now 



fit of the 1 
t no adaitio 
nd ranging d 
be included 
tern would b 
te responsib 
the pilot, 
ose calculat 
estimates. 



ow cost of 
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in the 
e designed 
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The micro 
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A. COMPONENTS 

The components of the modified system are depicted in 
Fig 4. Their relationship to the overall system is given 
below . 
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Fig 4 SYSTEM COMPONENTS 
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Control Stick Sense Switch 



1 . 



An on/off switch which is activated 'on' when the 
pilot begins tracking the target. 



2. Portable iiicro Computer With Auxiliary Memory 

A device whose auxiliary storage contains a complete 
library of standard release parameters and error sensitivity 
factors for each weapon authorized for carriage and release 
from the particular aircraft. Each weapon would require 
approximately 128 bytes of storage. It is assumed that no 
more than 50 such weapons would apply to any particular 
aircraft. The ground based micro computer would be capable 
of retrieving blocks of such data and passing the data to 
the micro computer aboard the aircraft. It is also assumed 
that no more than three different weapons would be carried 
by a particular aircraft on any flight. 

3 • Air Data Computer 

A device currently installed on Naval aircraft which 
provides analog signals corresponding to true airspeed, 
calibrated airspeed, and altitude. These signals would be 
converted to digital format to be read by the micro 
computer . 



4 . Loft Bombi np Computer 

An attitude reference gyro of the type currently 
installed in the F-4 whose analog signals corresponding to 
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pitch and roll angle as well as rate of pitch and heading 
change would be converted to digitaJL format to be read by 
the micro computer. 



5. Cockpit Keyboard 

Any standard numerical keyboard, of the type used on 
pocket calculators, which has at least five function 
switches and light emitting diodes to confirm register 
contents. This device would be used by the pilot, or 
preferrably by a crewman, to input the mode ( air-to-air or 
air-to-ground) , and if in the air-to-ground mode to input 
the desired delivery, weapon to be released, target 
elevation, and aircraft gross weight. 

6. Micro Computer 

A micro computer with a word size of at least eight 
bits and speed at least that of the INTEL 8080 (2.5 micro 
second instruction cycle) . 

7. Gunsigh t 

Either a standard fixed gunsight modified by 
replacing the fixed reticle light source with an array of 
addressable light emitting diodes, or a standard 'heads up' 
gunsight . 



8* Ai rb orne Mi ss ile Contro l System 

One of the several models of air-to-air radars 
currently installed on the F-4 which provides analog signals 
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corresponding to the present range to the target and the 
rate of closure. These analog signals would be converted to 
digital format to be read by the micro computer. 

9 • AngJ.e of Attack Sensor 

A device which is currently installed on Navy jet 
aircraft which provides an analog signal corresponding to 
the angle between the flight path of the aircraft and the 
FUSELAGE REFERENCE LINE. This signal would be converted to 
digital format for use in air-to-air calculations. 

B. DESCRIPTION 0? USE 

Fig 5 shows the overall timing sequence for the 
functions of the system. 
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1 . 



Air to Ground Mode 



Prior to flight, the aircraft is loaded with some 
authorized ordnance configuration. After the aircraft is 
operating on internal electrical power, a portable micro 
computer loads the appropriate tables into the memory of the 
aircraft's computer. As the aircraft approaches the target, 
the pilot/crewman keys a standard delivery and weapon code 
into the micro computer along with the target elevation and 
aircraft gross weight. As the pilot commences his roll-in, 
the micro computer samples calibrated airspeed and pitch 
angle and calculates the angle of attack of the fuselage 
reference line for a steady state dive at that pitch angle 
and calibrated airspeed and displays the angle of attack in 
mils on the gunsight. The pilot uses this information to 
establish the aircraft in a dive path which intersects the 
earth at some predetermined distance beyond the target ( aim 
off point) . When the aircraft is established in the proper 
dive, the pilot actuates the control stick sense switch and 
the computer begins to sample true airspeed^ calibrated 
airspeed, altitude, and pitch angle. The pilot makes no 
further throttle adjustments. The computer calculates dive 
angle and estimates what the release parameters { airspeed 
and dive angle) will be when the aircraft reaches the 
prescribed release height above the earth. This sampling 
and calculation takes about two seconds. The computer then 
calculates the effect of any errors in the estimated 
airspeed and dive angle and corrects for the impact error by 
addressing a new diode in the gunsight or sending a 
correction signal to the heads up gunsight. The new sight 
setting would be calculated using the error sensitivity 
factors which had been read prior to flight. The pilot uses 
the corrected sight setting as the new aiming reference and 
places the revised sight pipper over the target at release 
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altitude arid manually releases the weapon. After dive 
recovery, the computer is reset and the process is repeated 
as desired. 

2. Air to Air Mode 



For air to air use, there would be no need to load 
tables into the aircraft computer prior to flight. Instead, 
the program and appropriate tables would remain in the 
computer memory. 

In flight the pilot /crewman would select the 
air-to-air mode. The computer would then begin to sample 
true airspeed, altitude, present range to the target, 
closure rate, roll angle and pitch and azimuth rates of 
change . 

At this time the 35 mil setting would be illuminated 
on the gunsight. When the pilot had established a suitable 
tracking position with the sight held steady on the target, 
he would activate the control stick sense switch. 

The computer would then calculate the proper lead 
angle for the most recently sampled conditions and display 
this lead angle on the gunsight. The pilot would adjust to 
the new aim angle as smoothly and quickly as possible and 
fire the guns. Tracers would still be mixed with other 
ammunition so that the pilot could adjust his fire as 
necessary. The time required to perform the calculations 
would be about .2 seconds. 
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IV. 



ANALYSIS AND TEST 



A. AIR TC GROUND 



1 . General C on siderations 



If it is assumed that density, temperature, and the 
speed of sound are continuous, differentiable functions of 
altitude, then there exist unique solutions of equations (7) 
and (8) for each set of initial conditions. Moreover, it is 
intuitively clear that if the solutions were written in the 
component form 



0 0 0 0 

x(t) =F (v , v , x , y ,t) (9) 

x y 

and 



0 0 0 0 

Y (t) =H (v , v ,x , y ,t) , (10) 

x y 



then the position of the weapon at time t could be 
considered to be a continuous function of the initial 
conditions. That is, for any neighborhood, N, of a point 
which represents the position of the weapon at 



(x(t) , y (t) ) 
time, t. 



after release (with initial 



conditions 
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0 0 0 0 

(v , v ,x , y )), there exists a neighborhood, U, of the 
x y 

initial conditions such that any identical weapon released 
with initial conditions which were in U would, at time, t, 
be in N. Furthermore, it aav be assumed that first partial 
derivatives of ? and H exist within U so that 



dX = F * 
v 

x 

and 




x 



, 0 

* dV + 
v v 

y 



o 

F * dX + ? 

X v 




( 11 ) 



dY = H 



* 




+ 



H 



* 




+ H 



. 0 

* dX + 



H 



v 



X V 



X 



x y 

Where the variables and partials on 
changes in the initial conditions and t 
represent position at some time, t, which 



y 



the right 
hose on 
is held 



* c.'P. ( 12 ) 

represent 
the left 
constant . 



Since the proposed system incorporates no additional 
sensor devices, it rust perform essentially the sane type of 
calculations which the pilot of a fixed sight aircraft now 
performs. It samples airspeed, altitude, and pitch attitude 
and predicts what the value of true airspeed and dive angle 
will be at release altitude. These values are then compared 
with, the standard release parameters for the particular 
delivery being used. Any difference between the predicted 
value and the desired value is converted to a gunsight 
correction which must be applied to the sight to counteract 
the affect of the release airspeed cr dive angle error. The 
sum of the correction for airspeed and the correction for 
dive angle is then applied to the gunsight and the pilot 
flies the aircraft to place the new sight image over the 
target at release altitude. 

Since the proposed solution relies on the error 
sensitivity tables, and since these tables are based on the 
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partial derivative of F with respect to several of the 

components ot the initial condition ( for example, an error 

in true airspeed changes v ,v ,and x ) # and since the 

x y 

solution assumes that effects of changes in airspeed, dive 

angle, and aim angle are linear and independent, an 

analytical proof of the validity of these assumptions would 
require an explicit representation of F and H. 

Unfortunately, there is no known closed form solution 
available to cover all the possible sets of initial 
conditions. Therefore, a series of computer experiments 
were performed to determine whether the method of error 
sensitivity analysis was suitable. 



2 ♦ Te st s of the Proposed Solution 



a. Static Accuracy 

The first test of the proposed solution was 
designed to measure the effectiveness of the error analysis 
technique. To separate the error analysis technique from 
any errors which might be induced by airspeed prediction or 
dive angle calculation, it was assumed that these functions 
were performed exactly by the micro computer. To measure 
the ability of the system to improve accuracy, an arbitrary 
value of plus 200 feet or minus 200 feet was assigned to the 
uncorrected impact error, DX in figure 6. The original dive 
angle was assigned a value which was five degrees steep or 
five degrees shallow. It was assumed that the correction 
which was applied to the gunsight would cause the bomb to 
hit the target if the pilot could establish a parallel dive 
path whose displacement from the original dive path was -DX 
feet ( measured along the X axis). It was assumed that, 
instead of performing such a maneuver, the pilot used the 
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new sight image to moke a single change in his dive path. 
This change occurred three seconds prior to reaching release 
altitude and moved the aim off point -DX feet. This 
maneuver causes a small change in the dive angle for which 
the gunsight was not compensated. At release altitude, the 
effect of this change in dive angle was computed using the 
error sensitivity tables for the Mk-76 Practice Bomb (DX ' 
figure 6) . 



A FORTRAN program was executed to determine the 
values of DX', based on the previous assumptions. The 
program is listed in Appendix C ana results are in Table 1. 
In general the release error sensitivity method improved 
accuracy by a factor of eight, and, in the worst case, by a 
factor of five. 
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CORRECTION ALTITUDE 





DIVE AFTER CORRECTING 



. . UNCORRECTED BOMB 

CORRECTED BOMB 

Fig 6 SIGHT ANGLE CORRECTION IN A DIVE 
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TABLE 1 SIMULATION OF DTHET AND DX • 



b. Real Time Accuracy 



An important part of the proposed solution was 
the calculation of the angle of attack of the fuselage 
reference line. The pilot would use this information to 
establish the aircraft in a straight path which intersected 
the ground at a precomputed distance beyond the target, 
called the aim off point, reference/ - 3 The computer 
would also use this angle of attack information to compute 
the actual dive angle of the aircraft. 

Although all Naval carrier based jet aircraft 
are equipped with angle of attack sensors, these devices are 
designed primarily to provide accurate information in the 
landing phase of flight and are inherently inaccurate during 
high speed low angle of attack situations, ref erence/ 4 J. 

Therefore, in order to have a more accurate 
value for the angle of attack, figure 2-56 of reference/" 2 J 
was used as a source, - Fig 7. 
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DATA BASIS: FLIGHT TEST 



TAS-KNOTS (1 UNIT = 10 KNOTS) 



CAS-KNOTS (UNIT = 10 KNOTS) 




FUSELAGE ANGLE OF ATTACK - MILS 
(1 UNIT - 1 MIL) 



CAS-KNOTS (1 UNIT = 10 KNOTS) 



Fig 7 ANGLE OF ATTACK CHART 
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Figure 7 is the chart which a pilot would use to 
calculate the angle of attack for any particular air to 
ground delivery condition. The chart may be used with 
either of two sets of input keys. The first set is dive 
angle, temperature, true airspeed, and gross weight. The 
second set consists of calibrated airspeed, gross weight, 
and dive angle. Since the central air data computer 
provides calibrated airspeed information, the second set was 
used. To calculate manually the angle of attack for a given 
calibrated airspeed, gross weight, and dive angle the 
following procedure is used: 

(1) Find the desired calibrated airspeed 
(cas) on the right half of the bottom of the chart. 

(2) Move vertically to the curve labeled with 
the appropriate weight. 

(3) Proceed horizontally to the left until 
intersecting the straight line labeled by the desired dive 
angle . 



(4) Proceed directly down the chart and read 
the desired value of the angle of attack in mils from the 
bottom of the chart. 

To reproduce the information contained in figure 
7, it was necessary to represent the weight curves and 
straight lines in an algebraic form. The form chosen was 
the polynomial. Since all air-to-ground charts are based on 
the assumption that the aircraft gross weight is 42,000 
pounds with scalar adjustments for different weights, the 
same scheme was adopted for the micro computer system. The 
42,000 pound curve was approximated with a second degree 
interpolating polynomial, F (cas) . The straight Lines were 
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considered to be of the form 



y=m * x + b 

where ' y' represents ordinate values, ' x 1 represents 
abscissa values, ’m 1 is the slope of the line, and 'b' is 
the ' y' intercept value. Substituting ALPHA, the angle of 
attack, for ' x 1 , P (cas) for ' y', m (THETA) for m, and 

B (THETA) for 'b', the following equation was evaluated: 

ALPHA= (F (cas) -B (THETA) ) /m (THETA) . 

F was a quadratic interpolation of the 42,000 pound curve. 
The argument of F, cas, was taken as the difference between 
calibrated airspeed and 295. This linear translation 
allowed single byte raultiplicat ion to be used. THETA was 
the dive angle. B was a tabled, first divided difference 
function cf the zero mil intercept values of the straight 
lines in the lower left quadrant of figure 7. M, the slope 
of these straight lines, was found to approximate closely 
the value of '-SECANT (THETA) ' . 

Therefore, in order to speed the calculat ions, 
the equation which was actually solved was 

mil = cos (THETA) * (B (THETA) -F (cas) ) . (13) 

The COSINE calculation was based on a divided difference 
table which was arranged in four degree increments. 

An ALGOL program was executed to simulate an 
INTEL 8080 fixed point arithmetic evaluation of the 
proceeding equation. A listing of the program is contained 
in Appendix C. The results of the ALGOL simulation appear 
in table 2. A ELM version of the routine was timed by the 
INTERP8 as requiring 4798 micro seconds. 



37 



o 














CV? 


VO 


A 


A 


A 






A 


H 


i — 1 


i — 1 


▼H 


On 


CO 



o 














o 


A 


A 


CA 


CV? 






A 


tH 


i — i 


■» — i 


■H 


CO 


Cv. 



480 




A 

H 


A 


A 

i — 1 


12 


CO 


Cv. 


460 




\o 

tH 


A 

H 


13 


13 


ON 


oo 


o 

5 




CO 

i — 4 


Cv. 

■H 


A 

H 


t — ( 


o 

vH 


ON 


o 


cv 




H 


O 


o- 


VO 


■» — i 


O 






CV 


CV 


H 


H 


H 


■H 


o 


o 


_ 


VO 


A 


H 


o 


Ct 


A 






CV 


cv 


CV 


cv 


H 


tH 


o 


CO 




H 


o 


VO 




Cv- 


A 


CA 




A 


A 


CV 


CV 


i — 1 


H 



60 


NO 


A 


O 


ON 


O 


CO 


CA 


A 


A 


A 


cv? 


CV? 


■H 



o 

-4t 


A 


CV 


VO 


-3" 


-3' 


CN? 


A 


-3- 


-3- 


A 


A 


cv 


CV 



o 

cv? 


r— t 


ON 


A 


O 


ON 


VO 


A 


A 




-3- 


• 3 - 


cv 


CV 


O 

O 


i — 1 


ON 


T-l 


CO 


-3- 


H 


A 


vO 


A 


A 




A 


A 



n 

w 

CV 

a. 

co 


CO 


w 














cV 


ci 


M p 


O 


O 


o 


O 


o 


O 


M 

-u 




Q CV 
Q 


tH 


CV? 


(A 


-3“ 


A 


VO 



38 



T«i3LE 2. COMPJi'ED FUSELAGE KEFEttENCE LINE ANGLE OF ^TliiCK 



c. 



Forecasting Accuracy 



A third topic for investigation is the ability 
of the micro computer to forecast accurately the value of 
the release true airspeed early in the dive. The basic 
forecasting scheme is to fit a parabola through a series of 
samples of airspeed and altitude . It is assumed that both 
values can be read to the equivalent of four decimal places 
of precision. In fitting the parabola, altitude is treated 
as the independent variable and the release true airspeed is 
assumed to be 



2 

airspeed = C + C *relalt+C *relalt (15) 

12 3 

where •relalt' is the calculated release altitude. 

To determine whether such a method was 
appropriate, it was necessary to first construct a model of 
the F-4 in a dive, and then write a computer program to 
simulate the model. -The model chosen was that of Newton’s 
second law of motion: 

F=m*a 

where force, F, was assumed to be 

Thrust + W*sin (THETA) -drag. 

Thrust and drag were opposing forces acting parallel to the 
flight path, THETA was constant, and thrust, mass, and 
weight, W, were also constant. 

The resulting differential equations were 

A= (Thrust + K * SIN (THETA) - drag) /mass (16) 
and 
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V=A*t 



(17) 





Drag 


was assumed to 


be a function 


of 


and true 


airspeed. 


Figure 2-61 


of reference 


C 


selected 


as the 


source of data 


from which to 


obt 


of drag. 


This particular figure 


represents 


true 



vs. altitude at various dive angles for a particu 
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Runge-Kutta fourth order integration scheme. this was done 
in a FORTRAN subroutine, INTEG, (see Appendix E) . 

Using the model, three forecasting schemes were 
tested. The first scheme consisted of sampling airspeed and 
altitude values at one second intervals. A Newton divided 
difference interpolating polynomial was then generated and 
release true airspeed was predicted as 



tas 



release 



tas 



0 



+ (release alt-alt ) *F 

t (0 



+ (release 




( 18 ) 



alt-alt ) * (F 

t (0,1) 



F 



( 0 , 1 ) 



= (F -F ) / (X -x 
0 1 0 1 



/ 



F is analagously defined, 

0 , 2 ) 



F = (F -F ) / (x 

( 0 , 1 , 2 ) ( 0 , 1 ) ( 1 , 2 ) 0 




X represents the independent variable and F 
represents the dependent variable. 



The second method consisted of taking a total of 
twelve samples during a two second period. The first four 
samples were added and the average was taken as the true 
value. Similar averages were used to represent the two 
remaining inputs to a Newton divided difference second 
degree interpolating polynomial of the type described above. 



The third prediction scheme was to take 20 
samples during the same interval and to generate a least 
squares orthogonal polynomial to fit the data. The 
polynomial was also a parabola and it was generated and 
evaluated based on subroutines contained in reference /5 ? . 
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simulated. 



For the experiment, 1000 dive runs were 
Each run was based on the following assumptions: 

(1) Entry airspeed equals 345 knots. 

(2) Thrust equals 20,000 pounds. 

0 

(3) Dive angle equals 30 . 

(4) Entry altitude equals 6000 feet. 



(5) True airspeed and altitude would be read to 
four decimal digits of precision. 



(6) On any particular run, the altitude and 
airspeed error would remain constant. Eetween runs the 
error would _ be normally distributed with the mean airspeed 
error of zero and standard deviation of .5 knot. The 
altitude error would have a mean of zero and a standard 
deviation of 10 feet. 

The latter distributions are the result of 
assuming that the errors were normally distributed with the 
maximum allowable error equal to three times the standard 
deviation. Allowable instrument errors were obtained from 
reference/6 y. 

Normal distributions were simulated using the 
Naval Postgraduate School Computer Library package LLRANDOM. 
The test was repeated with a different seed to verify 
original results. See Appendicies A and E. 

d. System Performance 
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To determine whetner the use of the micro 
computer would significantly improve the delivery accuracy 
of the F-4 in the air to ground enviornaent, an experiment 
was conducted based on the following assumptions: 

(1) The previously described model represented 
an F-4 in a dive bombing run. 

(2) The ilk-8 1 Low Drag Bomb is a typical 
unguided air to ground weapon. 

(3) In the simulation, it would be assumed that 
the pilot, at release altitude, would have the sight aimed 
at the target. 

(4) All weapons would be released from a 
position 7 feet below and 18 feet aft of the gunsight. 

(5) All deliveries would be in no-wind 

conditions. 

(6) The path of the Mk-81 Low Drag Bomb could 
be described by the following set of equations which were 
obtained from reference 

X'' = - E * X ' (19) 



and 



y • ' = -E*Y' -g 

where X represents down range 
X', and Y' are the horizontal 
velocity vector. The value 
determined as follows: 



( 20 ) 

travel, Y represents altitude, 
and vertical components of the 
of the other variables was 
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E = (c * p * V * K ) /144 
a d 

where 



c is the reciprocal ballistic coefficient 
2 

(diameter / weight) , 



- (0.00003153 * Y) 

p is the air density 0.0751366 * e 



2 

V is the true airspeed in ft/sec . , 
a 



and 



K is the drag coefficient. 



The experiment consisted of simulating 1000 dive 
runs for each. of the following standard deliveries: 

0 



20 


dive. 


2000 


ft. release. 


450 


knots 


; true airspeed 


0 

30 


dive. 


3000 


ft. release, 


450 


knot 


true 


airspeed 


0 

30 


dive. 


5000 


ft.- release. 


500 


knot 


true 


airspeed 


0 

45 


dive. 


4000 


ft. release. 


450 


knot 


true 


airspeed 


0 

45 


dive. 


6000 


ft. release. 


500 


knot 


true 


airspeed 


0 

60 


dive. 


7000 


ft. release. 


500 


knot 


true 


airspeed. 



On each run the sampling method previously 
described as one second interval sampling was used. The 
predicted airspeed was calculated by equation (18). The 
entry conditions were considered to be normally distributed 



44 



about a correct set of (airspeed, altitude, thrust, and dive 
angle) which was calculated to achieve the desired delivery 
parameters. The standard deviations of the distributions of 
the entry conditions were 10 knots for true airspeed, 200 

0 

feet for altitude, 2.5 for dive angle, and 750 pounds for 
thrust. 



Normal distribu 
Naval Postgraduate School s 
deviations were assigned 
pseudorandom variables, x 
LLRANDOM package. The va 
uniform (0,1) distribution 
(0,1) distribution. The 
calculated as 

z= _DMEAN + y * STD 
desired mean, STD 
COIN is a statement 
x<.5 and a -1 other w 



tions were simulated using the 
ubroutine LLRANDOM. The standard 
with a data statement. Two 
and y, were obtained from the 
riable x was obtained from a 
and y was obtained from a normal 
desired value, z, was then 



V * COIN (x) where DMEAN is the 

V is the standard deviation, and 
function which returns a +1 if 
ise. 



After the release conditions had been predicted 
by the subroutine SAMP, two conditions were simulated. 
First, that the aircraft had no micro computer and that the 
pilot continued in a steady dive to release altitude. 
Second, that the pilot had a micro computer, that the micro 
computer supplied his gunsiaht with a corrected sight 
setting equal to the value calculated by SAMP , and that the 
pilot adjusted his aircraft's dive path as soon as this 
informa txon was available , so that at release altitude, the 
new sight setting was 'on' the target. 

For each condition, a simulated bomb was dropped 
at release altitude. A miss for the bomb was calculated as 
the distance between the point where the sight angle 
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intersected the earth and the down range travel, x, which 
was numerically computed ( based on the equations (19) and 
(20)). In addition, in the case of the modified flight 
path, a correct sight setting was also calculated. Finally, 
the ratic of the miss distance with the micro computer to 
the miss distance without the micro computer was calculated. 
If the micro computer gave a greater miss than the basic 
aircraft, the ratio was set to plus or minus 1.1, depending 
on whether the errors were both on the same side of the 
target. 



Arrays were built with these 1000 values and 
histograms were printed showing the distribution of the 
various values. See appendix B. The FORTRAN routine which 
simulated the experiment is in appendix F. An additional 
table was constructed from the histograms which compares the 
median miss distance under the assumptions that the micro 
computer tas or was not used. See table 3. 
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Micro Computer Suitability 



e . 



The previously tested errqr sensitivity 
algorithm was converted from FORTRAN to PLM for test on the 
INTEL 8080- The program was compiled on the Naval 
Postgraduate School IBM 360/65 and interpreted using the 
INTERP8 subroutine. Timing data was obtained and found to 
be .75 seconds for the entire program. See Appendix G. The 
program was then executed on the INTEL 8080 and final output 
(sight setting) was found to be in agreement with sample 
results which had been obtained from executing the FORTRAN 
version on an SDS 9300 computer. 



B. AIR TO AIR 



1 . General 

In air-to-air gunnery considerations, it must be 
assumed that the position of the target is some function of 
time, 

(x,y,z)= F (a (t) ,b (t) ,c (t) ) . ' N 

Since, in general, no assumptions can be made to 
restrict this motion relative to the (x,y,z) coordinate 
system of air to ground considerations, a new coordinate 
system is introduced. The origin of the new system, 
(x',y',z s ) is the muzzle of the fighter aircraft's cannon. 
The x' axis is always oriented so that a bullet's initial 
velocity vector will be aligned with the positive x' axis. 
The y' axis is perpendicular to the x' axis and forms, with 



the x* axis, a plane v/hich contains the wing roots of the 
fighter. The z' axis is . or thogo nal to the (x',y') plane 
with the positive direction measured in the same direction 
as the vertical tail of the fighter. 

It should be noted that the x' axis is not always 
aligned with the gun barrel boresight line. This is caused 
by the bullet's velocity vector being the vector sum of the 
basic muzzle velocity vector and the aircraft velocity 
vector. When the aircraft is flying at an angle of attack 
greater than zero these two vectors are not coincident. 
Therefore, the x' axis is shifted away from the gun 
boresight line toward the relative wind. 

The set of points from which a fighter can 
effectively employ guns against an airborne target lies 
within a 30 degree cone emanating from the tail of the 
target. The minimum firing range within this cone is 
approximately 500 feet and the maximum range is 
approximately 2000 feet. 

The dimensions of this firing region effectively 
restrict the problem so that, for a bullet to hit the 
target, the time of flight is less than two seconds. This 
short time of flight, together with the basic muzzle 
velocity of the standard aircraft cannon (3300 fps) , allow 
the computations of the proper lead angle to be performed in 
four distinct steps. Reference /$ J. The steps are as 
follows: 



(1) Estimate the future angular displacement of 
the target from the present line of sight. This 
displacement is assumed to be caused by the target's present 
velocity vector acting as a constant during the time of 
flight of the bullet. This angle is called the lead for 
target motion. 



0 9 



(2) Calculate the effect of gravity drop. 

(3) Calculate the effect of the fighter's angle 

of attack. 



between the 



(4) Calculate the 
gun and the gunsj.ght 



effect of 
(parallax) . 



the 



distance 



2 . Lea d for Target Motion 



The calculation of the lead for target motion in a 
plane parallel to the (y',z') plane is based on the 
following analysis of the situation: 

(1) If the fighter is tracking the target, 
that is, the fighter is in the firing cone with the 35 mil 
pipper held steady on the target., the relative positions of 
the two aircraft are 'close' to the positions needed for 
successful gun firing - . 

(2) Therefore, assume that the sight line and 
the x' axis are coincident. 

(3) At this time, the fighter's velocity minus 
the rate of closure represents the target's velocity 
component parallel to the x' axis. 

(4) The rate at which the fighter is turning, 
in radians/second, multiplied by the range is a close 
approximation to the target's velocity component parallel to 
the (y • , z ') plane. 

(5) If there exists a time, t', such that a 
bullet fired at time t-0 will strike the target at time t'. 
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then the geometry of the flight, of the bullet in a plane 
containing the x' axis and the component of the velocity 
vector of the target which is in the plane parallel to the 
(y« ,z») plane, can be approximated by Fig 8. 
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Fig 8 AIR TO AIR GUNNERY DIAGRAM 
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(6) Therefore, the lead for target motion 
orthogonal to the x* axis can be approximated by the angle, 
psi, (see ?ig 8) . 

The calculation of the angle in mils uses the fact 

that 



lim sin (x) /x =1. This implies that, for small angles, 
x-0 



ARCSIN (x/y) = x/y = AECTAli (x/y) 



where = means "is approximately equal to". 

In fact, for angles smaller than 200 mils, the maximum error 
of these approximations is 4 mils. Since the set of 
possible answers lies within the range 0*200 mils, the 
following equation was solved to determine psi. 

psi= VO/VX'. 

VO, the component of the target's velocity vector which was 
orthogonal to the X* axis, was calculated as the product of 
the fighter's rate of turn multiplied by the range to the 
target, multiplied by 1000. VX', the average velocity of 
the bullet, w<is calculated as the range to impact divided by 
the time of flight. 

Since the calculation, at this stage, ignores 

gravity, and since the bullet otherwise satisfies the 
assumptions of the air to ground weapon , provided that the 
x' axis is parallel to the surface of the earth, equation 
(19) can be modified to give 

A = - E * V (21) . 

x ' x ' 
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Since the velocity of the bullet during its time of 

flight will be several times larger than the speed of sound, 

and since the shape of the K vs. Mach curve is almost flat 

d 

above 1.6 times the speed of sound, equation (21) can be 



further simplified to 



2 

A = K *V (22) 

x ' x ' 

or equivalently, since V > 0, 

x ' 



2 

A / V = -p * K. 

X • X ' 

Since the effect of gravity is ignored, the density, p, will 
remain constant so that (22) becomes 

2 

A / V = -K , a constant. (23) 

x ’ x 1 2 

Integrating both sides with respect to t gives 

-1/V = -K * t + C . (24) 

x' 2 1 

Solving (24) for the initial conditions of the bullet as it 
leaves the gun gives 

1/V = (K * V * t +1 )/V (25) where V was taken 

x 1 2 m m m 

to be 3300 * V , and V was the fighter's velocity in feet 

f f 

per second. Expressing (25) in an equivalent form yields 

V = V / (K * V * t + 1) . (26) 

x ' m2 m 

Integrating equation (26) again with respect to t 
gives the equation for distance traveled by the bullet 
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x' = In (K * V * t + 
2 n 

Since x' at t=0 is zero, C 



C =0 
3 



1)/K 



2 



+ C . (27) 

3 



The time of flight may now be calculated as the 



time, t , when 
tof 



SANGE + (V 



f 



V ) * t 

c tof 




RANGE is the distance to the target at time t=0. V is the 



equation (28) is strictly greater than 1, a Newton's fixed 
point iteration function may be used to find the value of 

t . The iteration function is 
tof 

t = t - F(t )/F* (t ) . (29) 

n+1 n n n 

F(t) and F! (t) were defined as follows: 

F (t) = RANGE + (V - V ) * t -In (V * K * t + 1) /K (30) 

■F r* m2 2 



F ' (t) = ( V - V ) - V /(V * K * t + 1). (31) 

f c mm2 



The first approximation of t was used to find the 

tof 

values VX', and subsequently psi. To reduce the time 



c 



closing velocity at time t=0 



Since the argument of the natural logarithm in 



and 
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required to calculate t , t was always assigned the value 

2 1 

one . 



The Effect of Gravity Drop 

In calculating the lead for target motion, the 

effect of gravity uas ignored. A true trajectory, however, 

would always curve toward the earth as described by equation 

(20) . For a time of flight less than two seconds, equation 

<20) may be simplified by assuming that E=0. This allows a 

closed form solution for the distance travelled by the 

bullet below the no gravity flight path during time, t , 

tof 

seconds. The solution is 



2 

y = 32. 17 * t /2 (32) 

The error in (32) is less than 20 feet reference £ 8 J . 



The angle, in milliradians, subtended by the 

distance y , at a future range of x‘ feet is 
tof tof 



mils =y * 1000/x' 

gravity drop tof tof 

4 . Effect of AngJLe of At ta ck 

As the bullet leaves the barrel of the gun, it 
aligns with the relative wind. That is, the velocity vector 
of the bullet shifts from the gun boresight line to the x' 
axis. The amount of this shift, in mils, is calculated as 
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MILS 



= V * AC A/ ( V ♦ 3300) , 

AOA f f 

where AOA represents angle of attack, in mils, of the gun 

barrel boresight line. 



5 • Para ll ax Cor rection 

Since the gqn and the sight pivot are approximately 
seven feet apart in the z* direction, a small factor must be 
added to the sight setting. In the F-4, this setting is 
3.19 mils, reference /&/ . Fig 9 shows that the effect of 
adding a parallax correction is to have the zero sight line 
(ZSL) intersect the gun boresight line (TBL) at a known 
distance (harmonization range) in front of the fighter. 




Fig 9 PARALLAX CORRECTION 
6 . I tnp le mentation and Timing 

The calculations described above were coded into a 
PLM subroutine and tested on the INTERP8 simulation 
subroutine. The time of execution was found to be less than 
.2 seconds. In addition, a simulated sampling and data 
preparation routine was timed at approximately ,15 seconds. 
See appendix H. 



\ 
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V 



CONCLUSIONS 



Based on the experiments performed and the ti 
obtained, the micro computer is functionally 
providing an improvement in air to ground delivery 
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Althouqh no live tests were conducted 
either mode, such tests could be conducted 
low' cost by programming an A-6E or A-7E air- 
simulate the accuracy and timing of the 
computer system. This would elimina 
reconfigure the hardware of the test aircra 
the preliminary stages. 



on the system in 
at a relatively 
borne computer to 
proposed micro 



te 


the 


need to 


ft. 


at 


least in 



The cost of the micro computer and associated 10,006 
bytes of memory necessary to execute the air to ground and 
air to air programs is one thousand one hundred and fifty 
dollars based on currently available prices. Overall, the 
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micro computer r 
increase in the 



presents one way to help curb 
ost of aircraft procurement. 
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DISTRIBUTION OF ERRORS IN PREDICTED AIRSPEED USING ORTHOGONAL POLYNOMIAL 
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ERRORS IN PREDICTED AIRSPEED USING AVERAGED INPUT 
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ERRORS IN PREDICTED AIRSPEED USING ONE SECOND SAMPLES 
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SRRO B IN OOwN RANGE TRAVEL 
DESIREO OELIVERY 

JIVE ANGLE TRUE AIRSPEEO RELEASE ALTITUOE 

20. 450. 2000. 
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ERROR IN CCWN RANGE TRAVEL 



FREQUENCIES SAMPLE SIZE * 1000 
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RRCR IN UNCORRECTED DCrtN RANGE TRAVEL 



FREQUENCIES SAMPLE SIZE * 1000 
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FREQUENCIES SAMPLE SIZE = 1000 
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CGRkLCT EO ERRCR/UN CORRECT EC ERROR 
OESIRED DELIVERY 

LIVE ANGLE TRUE airspeed RELEASE ALTITUDE 
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ERROR IN CCWN RANGE TRAVEL 
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ERROR IN UNCORRECTEO OCwN RANGE TRAVEL 
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FREQUENCIES SAMPLE SIZE = 1000 
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CORRECTED ERROR/UNCORRECTED ERROR 



FREQUENCIES SAMPLE SIZE *= 1000 
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ERROR IN CCWN RANGE TRAVEL 



FREQUENCIES SAMPLE SIZE = 1000 
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ERROR IN UNCORRECTEO CCWN RANGE TRAVEL 



FRECUENCIES SAMPLE SIZE = 1000 
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DIVE ANGLE TRUE AIRSPEEO RELEASE ALTITUDE 



FREQUENCES SAMPLE S1AE = 1 uuu 
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CORRECTED ERROR/ UNCORR ECT ED ERROR 



APPENDIX C 



STATIC ACCURACY SIMULATION PROGRAM 
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I E ( THE T -3 0 . ) 60,80,100 
DR AG=FX(B( 4,1) ,3(4,21 ,8(4,3 ),TAS) 

RETURN 

D45 = FX( B( 3,1 ) ,B( 3,2) ,6(3,3) ,TAS) 

0) 0=EX( B( 4, 1 ) , 8(4, 21,6(4,3), TAS) 

UR AG= ( D30-D45 ) * ( THE T-4 5 . ) * ( - 1 . / 1 5 . ) +045 





o 










in 










o 










+ 














o 






• 




o 






in 


— 


o 


— * 


t/> 


*— >» i— < 


00 


4* 


oO 


<c 


UOUOX 


< 


oOOO— 


< 


ft- 


<<I • 


ft- 


<c 


ft- 


•* 


I — V— ♦ 


•* 


l-h- • 


r* 




- ► i 


— 


► ► 1 




<o 




CO 




(O 


•* 


(OrO-ft 


*• 


<0(0# 


•» 


CO 


0* 


00 


•* ►(-* 


<— 4 




COOJ' • 


O — 


—4 00 • 




oca 


^*-*0 


vOCQ 


— — O 


OcQ 


CO •* 


CD CD in 


•-H ► 


COCQO 


in ► 




► ► 1 




► ► | 




•*C\| 


— — 


Ooo 


— *■— ' *V— 


►00 


O ► 


OOOOLLI 


nJ- ► 


tMOOLU 


o ► 


•-HCO 


► *-x 


-- <(\J 


*> -X 


0- »*H 


r-H w 


COOsJl— 


•»«w 


<ooi — 


'~ 4w ' 


•*□0 




O CD 




f rp 


o ► 


CQ CD •ft 


OO ► 


cocn-ft 


o ► 


Q> *«-% 
r-4 


— O 


*-4‘— 

i~4 


o 


in — 

H H 


* r« 


»— 4 r— 1 ID 


•» 


r— 4 «— 4 vQ 


-— »* *k 


• rO 


•X ►O 


• oo 


► *o 


• r—4 


in'—* 


nN 1 


o— 


r-lOO 1 


O — 


CD 


^*—in 


in cd 




vOQD 




OOCQ^* 


1 — 


coaDin 


1 -* 




^ — Q^H-X 


— 


ft- X 


2:l^u.2:xx— 


o:ujilzxxw2iuu.z 
it ii , x ii ac 



■ r- ^ n * . ,, U. O- H U.-A- " _ 

53 1 ! || OHHOD II IIODHODQ 

r. UjN^inaro^LLCCUjNO^C^LUU.^U-!^ 
ZIOCOQO hQIOOO^^Q^ 



O o 

OO OJ 



oo 



o 

in 



o 

sOO 
•-H f** 



93 



APPENDIX F 



DELIVERY ACCURACY COMPARISON 
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CALL DIV( .DUM1, .DUM2, . F 1 0 i 

A3: 

/* MF10 := -F 1 0 */ 
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CALL SETSUP; 
L 00 P 2 : CAS- 100 ; 



THET = 30 ; 

K= MIL(CAS t THET) 
L00P3 : THET = 28; 

TAS= DI T A ; 

AL T= D I T A ( 1 ) ; 
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=MILS; 



EXEC PLM80CG, REGION. PASS 1= 1 50K , REGION. PASS2=150K 
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CALL MUL( X»Y,Z) ; 

RET UR NJ 

END MULT; 

COMPARE: PROCEDURE (XA-YA) BYTE; 

/* ROUTINE TO COMPARE TWO FLOATING POINT VARIABLES */ 
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